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Abstract. Recent developments based on relativistic hydrodynamic models in high 
energy heavy ion collisions are discussed. I focus especially on how hydrodynamics 
works at RHIC energies and how one can use the most of it in analyses of jet quenching 
and thermal electromagnetic radiations. I also comment on improvement of initial 
conditions and viscosity in hydrodynamic models. 



1. Introduction 

First data reported by the STAR Collaboration at RHIC pQ has a significant meaning 
that the observed large magnitude of elliptic flow for charged hadrons is consistent with 
hydrodynamic predictions [2 . This suggests that large pressure possibly in the partonic 
phase is built at the early stage (r ~ 0.6 fm/c) in Au+Au collisions at ^snn — 130 
and 200 GeV. This situation at RHIC is in contrast to that at lower energies such as 
AGS or SPS where hydrodynamics always overpredicts the data. Hadronic transport 
models are very good to describe experimental data at lower energies, while they fail to 
reproduce such large values of elliptic flow parameter at RHIC (see, e.g., Ref. |3J). So 
the importance of hydrodynamics is rising in heavy ion physics. In this review, I discuss 
first how well hydrodynamics can fit the data in Sec. 121 and next how hydrodynamics 
is related with other non-flow phenomena such as jet quenching and thermal photon 
emission in Sec. El In Sec. |H I argue several attempts to incorporate initial conditions 
from other models. Finally, I briefly mention about viscosity in Sec. 



2. Results from the hydrodynamic approach at RHIC 

After the first STAR data were published |T] , other groups at RHIC have also obtained 
the data concerning with flow phenomena [I]. To understand these experimental data, 
hydrodynamic analyses are also performed extensively j^j. Here I pick up several results 
on elliptic flow and emphasize on the momentum regions where hydrodynamics gives a 
good description. Elliptic flow is very sensitive to the degree of secondary interactions 
[Oj. The indicator is the second harmonic coefficient of azimuthal distributions [Zj 
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(For recent progress of higher harmonics, see Refs. (HUH]-) in ideal hydrodynamics, the 
mean free path of particles is assumed to be zero. So the hydrodynamic prediction of 
v 2 should be maximum among transport theories. The scaled elliptic flow, which is 
defined as v 2 divided by initial spatial eccentricity e = (y 2 — x 2 ) / (x 2 + y 2 ), becomes 
almost constant around 0.2 from hydrodynamic simulations fOj- Interestingly, the 
experimental data reach this hydrodynamic limit for the first time in central and semi- 
central collisions at RHIC energies [11J. The data of scaled v 2 at various collision energies 
monotonically increase with multiplicity per unit transverse area (1/ S)dN/dy and, 
eventually, comes to the hydrodynamic limit. The key quantity to achieve thermalisation 
is the number of particles in a unit volume. This can be seen also in the centrality 
dependence of v 2 . Deviation of the hydrodynamic result from experimental data starts 
from iVch/iVmax ~ 0.5 which corresponds to an impact parameter b ~ 5 fm [2j. The 
Pt dependence of v 2 contains rich physics. In low pt region, v 2 increases with pr 
almost linearly. On the other hand, the data points of v 2 saturate at high p T jT2] and 
deviate from hydrodynamic predictions. In the intermediate pt regions, the interplay 
between soft physics and hard physics is very important in understanding v 2 for identified 
hadrons. This will be discussed in the next section. Hydrodynamic predictions give good 
agreements with experimental data of v 2 including mass dependences in low pt region 
[IB*] . In low px region at midrapidity, hydrodynamics reproduces experimental data 
very well. While, in forward/backward rapidity regions where multiplicity becomes 
small, hydrodynamics overpredicts largely ^1] the experimental data observed by 
the PHOBOS Collaboration [15J. This suggests thermalisation is achieved only near 
midrapidity. The regions where hydrodynamics works well for charged hadrons are 
summarised as follows: | b \ £ 5 fm, j>r~ 1.5 GeV and | r] \ £ 2 jHJ. The scaled 
v 2 suggests that thermalisation is only partially achieved in small (1/ S)dN/dy. Heinz 
tried to discuss a possible mechanism of deviation between hydrodynamic calculations 
and experimental data in forward/backward rapidity regions or in peripheral collisions 
by introducing "thermalisation coefficient" [IB] . 

3. Information inside fluids 

Due to strong interaction among secondary particles, hadrons can be emitted only 
from freezeout hypersurface. So hadron spectra reflect the information simply about 
accumulation of the space-time evolution. Unlike the blast-wave model fitting, the 
space-time dependences of thermodynamic variables are obtained in hydrodynamic 
simulations. These informations are significantly helpful to understand what happens 
inside the reaction regions. Here I discuss two phenomena, i.e., jet quenching and 
electromagnetic radiations, which are directly related with information inside bulk 
matter. 
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3.1. Jet quenching 

Minijets produced in initial semihard collisions have to traverse the reaction region 
where bulk matter evolves. During traveling, these minijets lose their energies through 
interactions with the medium. So high px hadrons originated from fragmentation of 
minijets contain information about the medium [T~ • Gyulassy et al. first employed 
results from (2+l)-dimensional hydrodynamic simulations in a jet quenching analysis 
|18j . The first order term of an energy loss formula in the opacity expansion becomes 



Here kinematics of emitted gluons are neglected. A dimensionless parameter C includes 
strong running coupling constant and colour Casimir factors. One needs the space-time 
evolution of parton density p in quantitative analysis for the energy loss of a parton. 
Note that the parton density appeared in Eq. (J2J does not need to be thermalised because 
it simply comes from the inverse mean free path of an energetic parton. Nevertheless, 
the parton density in Eq. (j2J) should be a thermalised one at RHIC energies from the 
analyses of i> 2 discussed in the previous section. In Ref. [211], it is assumed that the 
parton density in an energy loss formula is a solution of hydrodynamic equations in full 
3D space which is compatible with low px data such as dN/dy and pt spectra 21J (the 
hydro+jet model). Systematic studies based on the hydro+jet model with Eq. (j2J are 
performed for p T spectra f22J, back-to-back correlation functions [2H], pseudorapidity 
dependence of nuclear modification factors |24j . and v% |22| 124"] in high px regions. It 
is not obvious where hydrodynamic px spectrum eventually turns into pQCD power 
law spectrum. Moreover, the transition point can depend on particle species. Interplay 
between radial flow and jet quenching is discussed in Ref. [22] • It is found that the 
transition point increases with hadron mass due to mass dependent effects of radial 
flow. The resultant vi divided by the number of constituent quarks is represented in 
Fig. 1 (left). This looks very similar to the results from recombination/coalescence 
models [221 • The nuclear modification factor Rcp for mesons is found to be smaller 
than for protons even though > m p |2E|. This may be the meson-baryon effect [22] 
rather than the particle mass effect resulting from radial flow. But there is a possibility 
that mesons do not participate in the hydrodynamic flow. This is suggested by a 
hybrid model in which hadronic afterburner is described by RQMD 27 . Theoretically, 
one needs to study hadronic afterburner by combining a hadronic cascade model with 
the hydro+jet approach. Experimentally, it is important to observe Vi for <fi mesons 
in very low px region where hydrodynamics is expected to work. These will reveal the 
production mechanism in the intermediate px region at RHIC. 

3.2. Electromagnetic radiation 

Photons and dileptons emitted from thermalised matter are signals of the QGP in heavy 
ion physics proposed many years ago [2H]- Here I focus on the discussion of thermal 




(2) 



Hydrodynamic models 



4 



c 

> 



0.2, 
0.15 : 
0.1 
0.05; 


-0.05E 



0.5 



— x 
.... K 




CGC+hydro+jet, b=5.5fm 



1 1.5 



2 2.5 3 
p T /n (GeV/c) 




7.5 -5 -2.5 2.5 5 7.5 



x[fm] 



Figure 1. (Left) V2 divided by the number of constituent quark from the hydro+jet 
model. (Right) Initial transverse energy density distribution of one specific event in 
a central Au+Au collision at ^snn = 200 GeV from NeXus. Courtesy of SPhcRIO 
collaboration. 



photon emissions. As emphasised by Moore [21], the yield of thermal photons is much 

more sensitive to temperature profile and space-time volume than the production rate 

(the invariant yield per unit space-time volume). In a hydrodynamic approach, the 

invariant yield of thermal photon is evaluated by accumulation of the production rate 

over all fluid elements above the kinetic freezeout temperature 

d TV ~ d¥} 

E ^ = E WE^iTizi), fifa), E = (3) 

Here E and p are, respectively, energy and momentum of an emitted photon measured 
in a local rest frame. Information from hydrodynamic simulations is taken through 
temperature T, chemical potential p (or fugacity A), four fluid velocity u^ 1 and four 
dimensional volume AVi for each fluid element at x,. A novel calculation of thermal 
photon yield is based on the hydrodynamics with rate equations for quarks and gluons 
[313 EI]. At collider energies, the QGP at the very early stage is supposed to be the 
"gluon plasma (GP)" since multiparticle production is dominated by small x gluons in 
the wave function of colliding nuclei at ultrarelativistic energies. Initial gluons achieve 
a thermalised state, but it may not be a chemically equilibrium one yet. Through the 
process gg — ► qq, the system of the GP goes towards to the QGP. Due to the smaller 
degree of freedom, the initial temperature in the GP is larger than in the QGP when one 
assumes a fixed initial energy density. The photon yield would be enhanced due to higher 
initial temperature. On the other hand, small values of fugacity A for quarks and gluons 
make the photon yield per unit space-time volume suppressed. The competition between 
these two effects is studied in terms of a hydrodynamic model with rate equations. It 
is found that unexpected cancellation happens between initial temperature effect and 
fugacity effect in the final spectrum . I note that it is worth studying thermal photon 
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yields from chemically non-equilibrium hadronic phase in order to conclude real excess 
of thermal photon yield from the QGP phase. Related with chemical non-equilibrium 
properties, it is interesting to see whether hydrodynamic results for elliptic flow will be 
changed at RHIC due to the insufficient chemical equilibrium. 



4. Improvement of initial conditions in hydrodynamic models 

All the hydrodynamic results mentioned in Sec.EJare based on ideal hydrodynamics with 
parametrised initial conditions which lead to the reproduction of the multiplicity and 
Pt spectra for hadrons in low p? region. Towards description of heavy ion collisions 
from initial colliding nuclei to final spectra in a unified way, it is much better to 
employ effective theories/models which are relevant in the very early stage of collisions. 
There are some attempts to incorporate results from event generators, e.g., HIJING 
[32] , URASiMA [33], VNI [33 or NeXus [3S], into initial conditions of hydrodynamic 
simulations. In other approaches, effective theories, e.g., string ropes/flux tubes |36j, 
the final state saturation [3Zj or the initial state saturation [38J, are employed. Here I 
pick up two approaches which are, to my mind, important for the RHIC physics. 



4-1. Fluctuation of initial conditions 

In conventional hydrodynamic simulations for non-central collisions, one parametrises 
initial transverse profile by using the number density of participants or binary collisions 
at a fixed impact parameter. This gives us a smooth profile for thermodynamic variables. 
However, if one picks up one event, the energy density distribution in the transverse 
plane has an extremely bumpy structure shown in Fig. 1 (right). In SPheRIO [3SJ, 
initial conditions for hydrodynamic simulations are taken from an event generator 
NeXus [33 as an event-by-event basis. The final spectra in this approach are calculated 
like conventional event generators: they perform hydrodynamic simulations, calculate 
particle spectra at each event, accumulate many events, and average particle spectra 
over simulated events. It is worth noting that the numerical cost of the hydrodynamic 
code SPheRIO is very cheap and that they are able to simulate many events even 
in full 3D space. It is found that the multiplicity from initial energy densities with 
fluctuations are always smaller than without fluctuations. One can easily show [10] 
the following relation between initial and final total entropy in a simple EoS case: 
('S') final ~ ('S') initial 1 — a ^e)^ wnere (' ' ') means the event average and a is a positive 
constant. One should keep this fact in mind in detailed analyses based on hydrodynamic 
models. Needless to say, a systematic study of the effect of fluctuation on t> 2 would be 
interesting since elliptic flow is sensitive to the initial geometry. 



4-2. Initial conditions from the colour glass condensate 

Jet quenching and large elliptic flow are two important findings measured at RHIC. The 
common key for these two phenomena is the dense partonic medium. What is a possible 
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origin of this dense matter in Au+Au collisions at RHIC? This can be traced to the 
initial parton density inside the colliding nuclei. It is well known that small x gluons 
are dominant for multi-particle production. In the ultra- relativistic limit, these gluons 
in colliding nuclei could form the Colour Glass Condensate (CGC) [JTJH2I- ^° mr tial 
conditions in relativistic heavy ion collisions can be described by melting the CGC [43J. 
Hirano and Nara employed a nuclear wave function discussed in Ref. [H] and calculate 
gluon number density at each transverse point through factorised formula. Regarding 
this gluon number density as a thermalised one, they evaluate energy density at each 
space point at initial time r . By throwing it into hydrodynamic simulations as an initial 
condition, they obtain pseudorapidity distribution in the usual hydrodynamic manner 
[35] and compare these results with the PHOBOS data [121 ■ The CGC initial conditions 
are found to work very well for describing the energy, centrality, and (pseudo) rapidity 
dependences of charged hadrons [3%] . 

5. Viscosity 

Finally, I briefly mention about the viscosity. There is no a priori reason why ideal 
hydrodynamics works so well at RHIC as discussed in Sec. El (For a perspective from a 
strong coupling gauge theory, see Refs. jlHl HZ]-) Viscous corrections to momentum 
distribution functions are considered in Refs. The first order correction 

to the momentum distribution function becomes 5f oc jifo(l + fo)p fl p u X^ u , where 
T s = §77/(e + p) is the sound attenuation length, f is the Bose distribution function, 
and X^ is the tensor part of the thermodynamic force. Relation between the inverse 
Reynolds number jSU] and the attenuation length becomes R^ 1 ~ T s /r in the Bjorken 
flow case |51j . R -1 is found to be very small from the blast-wave fitting with the 
above correction term |49j . This suggests that the hadronic fluid in heavy ion collisions 
is nearly perfect one. It is known that naive relativistic extension of Navier-Stokes 
equations breaks down due to infinite signal velocity since the equations are parabolic 
ones (see, e.g., Ref. |52j). One can introduce relaxation terms to avoid this problem |S2]. 
The existence of relaxation terms are essential since they change the type of equations 
to hyperbolic ones in which signal velocity remains finite. More systematically, one 
takes account of the second order viscous terms l>4, from the prescription of extended 
thermodynamics (see, e.g., Ref. Dynamical studies of viscous fluids are mandatory 

for comprehensive understanding of the QGP. 

6. Summary and Discussion 

I reviewed the hydrodynamic results at RHIC and discussed some related topics such 
as jet quenching and electromagnetic radiations from a hydrodynamic point of view. I 
also discussed recent attempts for initial conditions of hydrodynamic simulations and 
for viscous effects. In addition to the photon emission discussed in Sec. El dileptons are 
also interesting to study chiral properties of hadrons by using hydrodynamics. In the 
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low invariant mass regions, dilepton spectrum is the best tool to see the spectral change 
of hadrons in hot/dense medium (see, e.g., Ref. |56j). On the other hand, in the high 
invariant mass region, it is interesting to see whether J/ijj really melts in the QGP phase 
[SZlEHlEni- Hydrodynamics provides the temperature profile and can be also useful for 
quantitative analyses of these phenomena. Some open questions are as follows: Usually, 
the initial time in hydrodynamic simulations is chosen as around r ~ 0.6-1.0 fm/c. 
How do we get such an early thermalisation time? Although gg — ► ggg process seems 
to play an important role in thermalisation, the resultant thermalisation time is still 
a few fm/c [60J. Note that a typical life time of the QGP phase from hydrodynamic 
simulations is around 3 fm/c in central collisions at RHIC. In this review, I did not go 
into details about the HBT puzzle [SI] . The important point is to find a solution of the 
HBT puzzle which is compatible with other observables such as pt spectra and v<i- 
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